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TECHNICAL MEMORANDUM 


TORQUES ON THE GYRO IN THE GYRO 
RELATIVITY EXPERIMENT 

I. INTRODUCTION 


The purpose of the Gyro Relativity Experiment is to determine if the orientation 
of inertial frames is affected by the motion of nearby matter. General relativity 
predicts that local inertial frames will process due to the relative motion of the Earth 
as a whole (the geodetic precession) and due to the rotation of the Earth (the motional 
precession). The magnitude of these two effects in a low-Earth orbit are about 7 and 
0.05 arc sec/yr, respectively. It has been shown [1] that all viable metric theories 
wliich fall into the PPN formalism predict the same geodetic and motional precession 
as general relativity to about the 1 percent level. This report addresses the question 
whether the Newtonian drifts on the gyro as presently conceived in the Gyro Rela- 
tivity Experiment can be reduced to a level such that the geodetic and motional pre- 
cessions of general relativity can be detected. 

Section II is an extension of earlier works on suspension torques [2], A sim- 
plified expression is derived for the torque on the gyro for the orientation which is 
planned for the experiment, if the harmonic content of the rotor shape is known, 
our expression allows a simple hand calculation of the torques. This is an improve- 
ment over the more complicated general formula for torque given in Reference 2. 

This simplified expression allows a computation of the roll-averaged torque for all the 
harmonics. Some numerical results are given for some typical state-of-the-art rotors. 
Also, the effect of scratches on the rotor coating is considered. 

The remainder of the report is an analysis of the other torques on the gyro 
which are present. This is a parallel analysis to that of Everitt [3]. The results 
tu’e not much different, but the analysis was done independently, and many of the 
details of the derivations are supplied. Section III gives an analysis of torques due 
to gas drag from first principles with all the steps in the calculation supplied. 

Section IV gives an independent analysis of torques due to rotor charging, based on 
a result of Reference 2. Section V gives a more general treatment of mass unbalance 
torques than in previous work. Section VI considers gravity- gradient torques in an 
inclined orbit and exhibits a new expression for suspension torques caused by gravity 
gradients in an inclined orbit after roll and orbit averaging. Section VII gives a 
general summary of previous work on magnetic torques. Section VIII considers 
(xismic ray impacts and corrects a slight error 'n Reference 3. Torques caused by 
brownian motion, jitter, photon bombardment, patch effect, and dirt particles are 
not discu.ssed here, since they are small, and we have nothing to add to Reference 3. 


II. ELECTRICAL TORQUES 


The torques on the gyro caused by the suspension voltages have been analyzed 
in detail in Reference 2. In this section, we give a simplified expression for the 
torques valid for tlie planned configuration of the spin axis relative to the electrodes 


* i < ^ i = 
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in the experiment. We also calculate the effects of roll averaging on the complete 
torque expression including all the odd harmonics. 

The planned orientation of the spin axis of the gyro in the experiment is midway 
between two adjacent electrodes. That is, the direction cosines of the spin axis in 
the notation of Reference 2 are 


% = ^ ’ ^o " ^o " • 


Under these conditions the complete primary torque expression from pages 26 and 27 
of Reference 2 can be greatly simplified. Defining as before the three integrals 
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where the terms are defined in the list of symbols and in Reference 2, we can, by 
manipulating the limits of integration in the integral, prove the following identities, 
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Substituting these identities into the torque equations on pages 26 and 27 of Refer- 
ence 2, wo obtain upon using the expressions for acceleration f. and preload h. 
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where 



and the gyro spin axis is midway between the y and z axes and perpendicular to the 
x-axis. This expression is an improvement over pre\dous expressions because the 
dependence on f. and h. is explicit and all the harmonics are included in the a, b. 

and c expressions. The roll-averaged torque can then be computed directly from 
equation (1). Equation (1) has been checked with our previous computer program 
described in Reference 2 and agreement has been obtained. For a given rotor, we 
need only compute a, b, and c, and we have the complete torque. 

If we expand r(6') in a fourier cosine series, the coefficients of the a^'s can 

be computed numerically (as in Table 1). These numbers have been checked analyt- 
ically for the second, third, fourth, and fifth harmonic and found to be accurate to 
six figures, verifying the accuracy of our numerical integration scheme. We notice 
that even harmonics only contribute to b, and odd harmonics only to a and c. Also, 
the coefficients do not decrease much as the order of the harmonics increases, veri- 
fying earlier work. So if we have the harmonic coefficients a^ for a given ball, we 

can use Table 1 to compute the complete torque. In practice, the a 's will alternate 

n 

in sign as will the X^, Y^, and Z^, so a lot of cancelation will take place. The 

result of the summation for a. b, and c will then not be significantly greater than 
any one term in the series except perhaps in pathological cases. So we can get a 
good idea of the magnitude of a, b, and c by just looking at the general size of each 
term in the series. This gives us u way to quickly estimate the drifts for any given 
ball. 


To illustrate this, consider the three rotor profiles in Reference 2, Figures 1, 
2. and 3. The values of a, b, and c for these rotors are given in Table 2; we see 
that the general level of the harmonics is reflected in a, b, and c. b is dominated 
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TABLE 1. COEFFICIENTS OF (6^ = 30 deg) 


Harmonic No. 

Y. 

1 

Harmonic No. 


Z. 

_1 

2 

-3.84765 

3 

-2. 94524 

1.76715 

4 

0. 96191 

5 

3.68156 

-1.47262 

6 

2.52502 

7 

-0.27918 

0.25771 

8 

-0. 81161 

9 

1. 30465 

0.82835 

10 

2. 78992 

11 

1.50071 

-0. 96181 

12 

-0.73694 

13 

-0. 97497 

0.20939 

14 

0. 68016 

15 

1. 45899 

0. 62342 

16 

0.02978 

17 

-1.69550 

-0.77133 

18 

-1.46040 

19 

0.11291 

0.17999 

20 

0. 57729 





20 





r(£0 = Y. \ 

cos n0' 




n=2 





20 





a = V X a 
n n 





n=2 





20 





b Y Y a 

Z- n n 

n-2 


Harmonic No. 

X. + 2 Z. 
1 1 




3 

0.58906 


20 


5 

0.73632 


c = V Z a 
n n 


7 

0.23624 


n=2 


9 

2. 96135 




11 

-0.42291 




13 

-0. 55619 




15 

2.70573 




17 

-3.23816 




19 

0.47289 




TABLE 2. ROTOR TORQUE COEFFICIENTS (yin.) 



a 

b 

c 

Case I 

0.018 

4.3 

-0.078 

Case 11 

0.207 

14.5 

-0.106 

Case III 

1.065 

10.1 

-0.337 


by the large second harnronic, while a and c are dependent on the general level of 
the odd harmonics. Case III shows that higher harmonics can be fairly significant. 
However, the random nature of both coefficients in the series for a, b, and c improves 
the situation. Table 3 shows the harmonic coefficients for our sample rotors from 
Reference 2. 


TABLE 3. HARMONIC COEFFICIENTS (yin.) 


Harmonic No. 

Case I 

Case II 

Case III 

2 

-1.0974 

-3.6691 

-2.5278 

3 

-0.0797 

-0.0958 

0.0217 

4 

0.0578 

0.1878 

-0.0921 

5 

-0.0420 

-0.0328 

0.2589 

6 

0.0656 

0.0400 

0.2019 

7 

-0.1091 

-0.0445 

-0.0820 

8 

0.0897 

-0.0704 

-0.0999 

9 

-0.0150 

0.0052 

0.0312 

10 

-0.0022 

0.0366 

-0.0263 

11 

-0.0550 

-0.0055 

0.0069 

12 

0.0158 

0.0302 

0.0756 

13 

-0.0215 

0.0023 

-0.0536 

14 

-0.0151 

-0.0004 

0.0150 

15 

-0.0122 

0.0014 

0.0489 

16 

0.0301 

0.0010 

-0.0481 

17 

-0.0031 

-0.0209 

0.0123 

18 

0.0112 

0.0124 

0.0203 

19 

-0. 0092 

• 0.0032 

-0.0136 

20 

0.0037 

0.0035 

0.0160 
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A series of computer runs was undertaken for some general rotor orientations 
with the program described in Reference 2. These computations indicate that the 
general level of torque is given to gxxtd approximation by equation (1) for orientations 
near to that assumed in calculating equation (1). This means that the torque is not 
very sensitive to misalignment of the rotor spin axis in the experiment. 


In the final gyro experiment, the spacecraft will be rolled about the spin axis 
of each of the gyros with a period of about 15 min. So we need to calculate the roll- 


averaged value of equation (1). 
obviously zero if h 


Since the roll-averaged value of the i^y'h^ terms is 
and h^ remain constant over the roll period, we need only con- 


sider the terms involving the T's. We first define a new set of axes such that the z’ 

axis is parallel to the spin axis and the y' axis is perpendicular to it and the x ~ x* 
axis. In this coordinate system we have 


fy' = y ‘fy - 


f ’ - ~ (f + f ) 


T ^ -^—'5 

X „ . 2 „ 

2 .! sin 


f 'f ' 

f ' a + b 

y h 


a'f ' 

y 


T ' = -L (T -T ) 

y ^5 " y z 


m 


sin ” 


„ (f c) c' f 
2 X X 


0 


where we have assumed h = h - h. 

y z 

We can now average over the roll by considering a transformation from the 
rotating spacecraft coordinate system X', Y', Z'. to an inerti ally- fixed coordinate 
system X". Y". Z" 


X" X cos , t - Y' sin ^ t 


Y" - X sin . t + Y' cos . t 
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Z" = V 


OF 




where ui is the rotation rate. 

cx)S u)t “ 

Ty„ = sin wt + Ty, 


Then we have 
sin u)t 
cos 0 )t 


substituting in the inverse transformation 

T^, = a' [- f^„ sin ut + f^,, cos ait] 

T , = c' (f cos u)t + f sin ojt] 
y X y 

and averaging, we find 



which relates the roll-averaged torques to the accelerations in the inert ially- fixed 
system. The roll-averaged drift will then be 


A 

2tt 


a+2c+bf^/v15 h 


r 

o 



V 


( 2 ) 


for (ij = 30 deg where v is the peripheral velocity of the ball, f^^ is the magnitude of 

the acceleration perpendicular to the spin axis, and f. - f , is the longitudinal accel- 

eration. Thus, if we know a, b, and c for a given rotor, we can calculate .i for any 
acceleration level. The values of H for the three rotors considered in Reference 2 
are given in Table 4 for several values of fn/h. Drifts scale as f. . So we see for 

^ -10 

the balls consider, we are at about the 0 1 mas/yr level for f = 10 g and corre- 

^ - 8 

spondingly higher values for larger f^^. For example, f^^ =10 g, a level which 

may be attainable witlwut drag-free control, gives drifts of the order 10 mas/yr. 
Notice that the roll averaging does not reduce the odd harmonic torques much as can 
bo seen from Tables 1 and 2, where there is not much cancelation in the a + 2c terms 
except for the third and fifth harmonic. 
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TABLE 4. ROLL-AVERAGED DRIFTS (MAS/YR) 
(fjj, = 10 g V = 2400 cm /sec) 


fj/h 

Case 1 

Case II 

Case III 

10-^ 

0.037 

0.001 

0.110 

10'^ 

0.029 

0.026 

0.128 

10'^ 

0.032 

0.282 

0.306 


For rotor misalignments and h ^ h , we will get imperfect roll averaging of the 

y * 

even harmonics, and this effect will be proportional to the misalignment angle. The 
exact magnitude of this effect for all the harmonics has not been calculated, but for 

- 4 

misalignment angles of ^ 10 rad, the effects are negligible. 

In the process of depositing the niobium film on the gyro rotor, it is possible 
that small scratches on the film may occur. To estimate the effect of these scratches 
on the torque, we consider the following situations. Suppose the ball has a line 
along a great circle in which the ball radius changed by Ar in a distance At. This 
would represent one side of a scratch. Then along that scratch, we have 


dr( i') 

di7 


' 1’ ( r ) 
o 


the solid angle covered by the scratch would be 


(At)(2-y^) 


o 


Hence the contributions to and can be estimated as 


dr(^) 

"df' 


sin 


min 


r 


sin 


nun 


. 4 r 
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where ®min is the minimum between the gyro spin axis and an electrode = 

15 deg). For a 100- yin. scratch, this could be a large effect, except for the follow- 
ing points. The other side of the scratch will contribute an equal and opposite 
torque if the scratch is symmetrical. Furthermore, r(6’) is the average shape of the 
rotor over a spin cycle so that, tinless the scratch happens to be perpendicular to 
the spin axis, the contribution to the average shape is 


2tt 


Ar’ = 


/ 


Ar(6,4i)d<}) O' Ar 


(It) 


so the effective Ar is reduced by a factor (At/r^) due to spin averaging of the shape. 
For At ^ Ar 100 yin., the contribution to I becomes 

. g 

I ^ 4Ar' ^ 4 ' 10 yin. 


which is considerably less than the contributions to I^, and 1^, due to large- 

scale ball asphericity. Thus, we conclude that sufficiently thin or shallow scratches 
are not a problem, except under very improbably circumstances. 


III. TORQUES DUE TO GAS DRAG 


The residual gas in the gyro housing will cause rundown torques on the rotor 
if the pressure is uniform and drifts of the gyro if the pressure is nonunifortn. We 
calculate these torques in this section. Since the mean free path of the gas molecules 
is many orders of magnitude huger than the dimensions of the rotor-housing gap and 
the pressure is very low, we can use the kinetic theory of gases for an ideal gas. 

We calculate the number of molecules impinging on the rotor per unit time and then 
assume that each molecule of mass m transfers a momentum mV to the ball, where V 
is the peripheral velocity of the ball at the point where the molecule hits the ball. 
Integrating over the ball then gives the total torque. 

Of those molecules which have an x component of momentum p^, (N/V)dAtp^/m) 

of them will strike per second on an area dA, perpendicular to the x axis, N/V being 
the number per unit volume and (p^ym)dA being the volume of the prism, in front of 

dA, which contains all the molecules that will strike dA in a second. The average 
number of molecules striking dA per second is tlx^n obtained by integrating over the 
Maxwell distribution. 
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2m kT 


N ha 
V dA 
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where k is Boltzmann^s constant and T is temperatiire. Multiplying by mV and dividing 
by dA, the momentum transfer per unit area per unit time B is then given by 


B = V P 




m 


where F is the pressure, and we have used the ideal gas law. This agrees with 
Reference 4. equation (1.113) except for a factor of order 1. This factor accounts 
for the fact that the total momentum mV may not all be transferred to the ball, 
depending on the character of the surface of the ball. Our equation for B is clearly 
a worst -case calculation. 

If we define a coordinate system by unit vectors x, y, and z, where z defines 
the spin axis of the ball, then the momentum transfer will be in a direction perpen- 
dicular to the spin axis anJ tangent to the ball, that is, in a direction p given by 


= ii sin : - ^ cos 4' 


where 4 is the azimuth angle from the x-axis. The total torque is then 


T - J (r X . )BdS 


(3) 


where v is the radius vector to the surface element dS. Defining 0 as the angle from 
the 7 , axis and using V - .. r^ sin 6, wc find, rtssuming P is uniform. 


T 

7 . 
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2'rkf 
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this produces an exponential rundown of the form 


0 ) 


w e 
o 


-t/T 


where 


T 


1 /27TkT 

5 P P / m 


in agreement with Reference 3, equation (7), where p and r are the density and 

- 9 ^ 

radius of the ball. For P = 10 torr and T = 1.6K, we find t = 300 yr so the ball 
speed will decrease by about 0.3 percent in a year. From equation (3), we find 


T 

X 


m 


27TkT 




I 


sin 9d6P [sin 6 cos 6 cos 


o 


4 >] 


and a similar expression for with cos <p replaced by sin If P is uniform, both 

these expressions vanish. However, if P varies in the housing, it can cause drift- 
producing torques. Letting P = P + AP and letting AP vary over an area A, the 
worst -case drift becomes 


, _ 15 / m A P A 

“p 16 7T / 27TkT 3 

p r 

o 


where we have replaced the bracketed quantity in the integrand by its maximum value. 


2 

1/2. Tliis is similar to equation (126) of Reference 3. Letting A = 1 cm , A? = 

10 torr, and T = 1.6K, we find 1 arc sec/yr. So we see that drifts caused 


by pressure variations can be quite large. However, they are averaged very well by 
spacecraft roll, since the pressure variations will roll with the spacecraft. Since the 

averaging effect goes as 1/# rolls and at a 15-min roll period there are about lO"* 
rolls in a year, the drift is reduced to about 0.1 milliarcsec/yr which is similar to 
the number given in Reference 3, page 571. Again, we are relying strongly on space- 
craft roll to reach the accuracy goal of the experiment. However, if the gyro housings 
are sealed during the experiment, the pressure uniformity may be considerably better 
than the 10 percent assumed here. The effects calculated here could be tested on the 
ground at higher pressures and extrapolated to the lower pressures, but they could 
not be checked at the lower pressures required for the experiment on the ground, 
since drifts from other sources are much larger. Reference 3 reports data in which 
the rundown was measured to be about 1 order of magnitude larger than that expected 
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from the theory presented here. This extra rundown may have been due to the 
suspension system which can cause large rundown at high g-levels. This effect would 
be greatly reduced with a low-g suspension system. 


IV. TORQUES DUE TO BALL CHARGING 


The rotor in the gyro experiment will acquire some charge due to the charged- 
particle environment in space and possibly electrical breakdown during the initial 
levitation. In this section, we calculate the torque due to the electrical charge and 
put limits on its size. 

From Reference 2, the torque on the gyro is given by the expression 


f - 


0^0 


2d 


o 






yB^) + iya^ 


o.Y„) 


+ F. (ap. 




(dr(eO/d9^) 

sin 0’ 


sin0d9d({) 


(4) 


where the terms are defined in the list of symbols and in Reference 2. For a charged 
rotor, we can then take all the V.'s to be given by 


V. 

1 



o r 


(5) 


where Q is the total charge on the ball. This follows from Gauss' law if we assume 
that the charge is distributed uniformly to first approximation. Then, if o is the 

2 

surface charge per unit area, the radial electric field is E = <jJ e - Q/4tt£ r , 

r o o o 

following the construction of a Gaussian surface as in elementary physics. Since the 
rotor- electrode gap is much smaller than the rotor radius, we can assume a parallel 
plate capacitor and V. = E^d^ from which equation (5) follows. This equation also 

follows for the situation of two spherical conductors with the inner one charged, 
which our situation here only approximates. 

We can now use equation (4) to evaluate the torque. If we neglect the ^d/d^ 

terms and consider the case 3^ - = 1/v?, we can see from the equations on pages 

26 and 27 of Reference 2 that "f = 0. This follows from the identities in Section II. 

If we are not precisely at this symmetry point (which is the planned configuration for 
the experiment), then it is not obvious from the general expression that T - 0, but 
expansion in terms of harmonies indicates that it is, since only differences in voltages 
occur, rius means that the torques will come from the Ad/d^ terms in equation (4). 

Since these terms are very complex to evaluate (Reference 2, Section VIII), we will 
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only estimate their order of magnitude. Taking, for example, the contribution from 
the z+ electrode, we find 


2 + 


31 ^ 

~2 ^ 


* ( Jo + Jq) + — (■^1> ~ 

X. /HT X 


32n‘e r ‘ L ^ ^ 

o o 


V 2 


where 


27T 


= / -n (^) 
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dr(6’)/d0' 


sin 


0’ 


0, 2tt 

Jj = r sin 6d0 y d(|, (^g-j 0 -grg-gi 


2tt 


Jg = r sin 0d9 

o o 


/ (-t) 


Ad\ dr(e’)/d0' 

sin 6 ' 


eos 0 ' = — [sin 0 sin (f- + cos 0 ] 


Similar terms appear for the other five electrodes. An estimate of the order of mag- 
nitude of these terms can be obtained by factoring out (Ad/d^) and taking only the 
second harmonic for r(0’). For example. 


0. 2ti 

Ad . r . r , , [sin 0 sin + cos 0] q 

j ^ ^4a« / sin 6d9 / d([) i — cos 0 

^ d 2 J J 

o o 


Sna. 


1 - cos 0 , 




For Oj = 30 deg, then just for the Jg term on the z+ electrode, we get 
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(0.117) 


Multiplying by a conservative factor of 3 for the and J 2 contribution and a factor 
of 6 for the six electrodes , we find 


X 


Q‘ 


4 F r IT 

o o 


(Sf) 


This is a very conservative number since a lot of cancelation will take place if the 
original expression for ^ were calculated exactly. The higher harmonics will contribute 
to the torque also, but vye know from Reference 2 that is the dominant contributor 
to these calculations. This gives a drift rate 



which is similar to that of Reference 3, equation (120) [equation (120) requires a 

factor of or in the denominator to be dimensionally correct]. Now, assuming a« = 
o ^ 

2 uin.. Ad/d^ ^ 1/100, V - 2400 cm/sec, we have as a limit on Q in order to obtain 
0. 3 milliarsec/yr drift 


Q = 1 • 10* electrons 


which is the same order of magnitude as in Reference 3, page 561. We know, how- 
ever, that the Ad/d terms are averaged by roll so that assuming the standard reduc- 
- 4 ° 

tion of 10 due to roll averaging, we get a limit of 
9 

Q = 4 '10 electrons 


which is again similar to Reference 3. 

Whether this level of charge can be maintained is an open question. Reference 
3 gives some data which indicate that for the DISCOS satellite, the charge level was 
kept acceptably low. Reference 3 gives a discussion of this and a method to measure 
the charge buildup. 
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If the center of mass is displaced from the center of support of the rotor, then 
there will be a torque of the form 

T = (Ar) S X f , 


where S is a unit vector in the direction of the gyro spin vector and f" is the external 
force. Ar is defined as 


^ rrp(r)dV 

Ar = S • <- 

/ p(?)dV 


where the volume integral of the density is taken over the actual mass distribution. 
The origin is taken as the center of si^port of the rotor. For a spherical coordinate 
system with 0 defined with respect to S, the integral can be written 


2ti 


mAr = 


r(9,(|)) 


/«/ sin ed0 J 


r'^dr 


p(r,6,(|)) cos 9 


where r(9,$) is the actual shape of the rotor and m its mass. We wish to calculate 
the effect of nonuniformities of p and r(6,4)) from this integral They will lead to 
worst -case drifts of the form 



We can expand p in spherical harmonics [5] 

0° i 

Pir.e,*) = Z L Y^^( 0 ,« . 

!l=0 m=- ?, 


If we now consider only density nonuniformity in a spherical ball, then, since 
Pj(cos0) = /4tt/3 Yj^q(0,4)) = cos 0 and using the orthogonality properties of 
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Or 


mAr = 


4TT 


where 




3 

j r dr PjQ<r) 


o 


2tt tt 


Pio(r) 


f J' P(r,9,4t) cos 6 


o 


Given the actual density profile, then this integral could be evaluated. Notice that 
only the lowest spherical harmonic coefficient contributes to the torque. If we now 
assume PjqCi*) is independent of r, then, for a given drift rate Q , this requires 


-''lO 8 V 

p ^5 / 3 r ““o 


6 V 1C'® 


for V = 2400 cm/sec, f = 10 g, and = 0.5 milliarcsec/yr. This is a condition 
on t Yo the higher If density variations were due only to Pj^(cos 6), then 


A, p(9=0) - „(0=rt = /I , 

but this is an overly conservative criteria and much larger peak -to -valley excursions 
in density can occur, provided they can be characterized by spherical harmonics 
other than 

If we had a bubble at (R^, 0^, 4)^) , then 


m'r = ■' m R cos 0 
o o 


and for the above parameters 


\m ^ o 1 

m R cos 0 
o o 


10 



or the bubble has radius 


r 0.2 mm 
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Now consider a nonspherical ball of uniform density and radius 


r(e,(})) = + Ar(6,(l)) . 

If we choose the origin at the center of the ball for Ar = 0, then the r^ term con- 
tributes nothing and 

2it 7T 

-A. = / sin 6d6 Ar(6,4)) cos 6 , 

o o 

where we are considering Ar as due to a coating of density p^. Expanding Ar in 
spherical harmonics and using the orthogonality relations again gives 


Ar(e,4,) = E E 

£=o m=-£ 


mAr = 


4tt 3 
3 ^c^o ^10 


where 


‘10 



2tt 7T 

d(^ J sin 9d0 Ar(6,4i) cos 
o o 


or for a given drift rate 


< I A V _r ^ ^ 

r 5 / 3 f p„ o 


-7 


for 4 where is the rotor density and the previous values of V, f, and 

If Ar(0,(t)) is known, we can compute a^^Q fr-orn the above integral and check the 
condition on a^^. If Ar(fi.(|)) were all due to Pj^(cos 0), then 
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Jii 


Ar(&=0) - Ar(9=TT) = /-a-^ 

TT 10 


but again this ir> over pessimistic. In particular, for ball-coating procedures in which 
selected, equally- spaced sites are coated, the coating shape would be largely due to 
the higher which do not contribute to Ar. So much larger density variations 

could be tolerated than are evident from the condition on a^Q. A more quantitative 
statement requires an actual coating profile. 

Finally, we consider a scratch of depth Ar’ and width d = r^AB at 9=6^. Then 


2 

mAr = 27ip r sin 0 cos 0 (Ar'd) 
CO o o 

which requires for previous values of 


O 

im 

r 

o 


4 
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P sin 0 cos 0 
c o o 


(1.9 


10 


For Ar' = d = 100 pin., this seems to be satisfied for P^,/Pj. 4. So from the point 

of view of mass unbalance, we can tolerate fairly large scratches on the rotor. 


VI. GRAVn : GRADIENTS IN AN INCLINED ORBIT 


In an inclined orbit, the spin axis of the gyros and spacecraft will not be in 
the orbit plane. This introduces extra terms in the gi’avity- gradient expressions not 
considered in Reference 2. Section IX. These terms are discussed in this section. 

If the orbit is in the y-z plane and the spacecraft and gyro spin axis make an 
angle with this plane, then using the expression for the gravity- gradient forces in 
Reference 2, Section IX. the gravity- gradient forces are seen to be 


-0 cos u,t cos 6 - D sin 0 
o z 

0 

(3 cos" ' t-1) p cos a t+3 cos sin P't(-p^ sin (jjt sin 0+p^cos B) 

2 

(3 cos d sin ;t) o cos ,t + (3sin ;d-l)(-p sin ;r,t sin B+p cos 9) , 
o o z 


where 

. ^ - displacement of gyro perpendicular lo spin axis of spacecraft 
, displacement of gyro parallel to spin axis of spacecraft 
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IS 

OK POOR QUALITY 

0 ) = spacecraft roll frequency 
= orbit angular frequency 
R - orbit radius . 


If we now average this over time in the manner of Reference 2, Section IX, we find 



Transforming this to a coordinate system with the z' axis coinciding with the space- 
craft spin axis, we find that the component of force perpendicular to the gyro spin 
axis is 


<f, > = 

tr 


3GM 


2R 


3 ^z 


p sin 0 cos 


and the force along the spin axis is 


. . GM / 3 2 ,, , 


( I ) 


_ g 

For 10 cm, 10 g» this will produce torques on the gyros of order 

10 mas/yr from Table 3 with our sample gyros. This is 2 orders of magnitude larger 
than gravity- gradient effects in a polar orbit, indicating that a polar orbit is a much 
better choice for the experiment. The only other possibilities would be to (1) mini- 
mize in the "four square" configuration or (2) use the known time dependence of 

P which can be calculated from the known orbit regression rate to separate the 
gravity- gradient drifts from the relativity drifts. Considering the effort involved in 
this experiment, neither of these alternatives seems reasonable. So, it would seem 
that these gravity- gradient considerations make a polar orbit imperative. 

Notice also that the product f^ ^tr^^ equation (2) can be fairly large and 

does not average to zero for the gravity gradients. For h 10 ^ g, this term is 
comparable to the other terms in equation (2), so it must be included in the analysis 
if modeling is attempted for the torques in an inclined orbit. The exact form of this 
term after orbital and roll averaging equation (1) in the manner of Sections I and IX, 
Reference 2, is 
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<T. 


-b 


16/J h 2tt sin 


m /gM 2\ 

. 2„ { p3 ) 

in yR J 


4 3 

(cos 0 ' sin 0 cos ^ + 3 sin 9 cos 0) 


<T 


b 


11 


4 h 2 


f- ^ ( 

TT sin^Oj^ yR** / \ 


2 / 27 2 

- I 3 - cos 0 I sin 0 cos 9 


sii 




2 3 2 2 3 

2 cos 6 ■ 2 ® ® ® ® 


) 


where means the component parallel to the orbit plane and means the component 5 

perpendicular to the orbit plane, p is the displacement of the ;?yro from the spin 
-62^ 

axis. For h 10 g, the p term leads to drifts of order 10 mas/yr for p 'v 10 cm. 

z z 

Thus the "in line" configuration produces unacceptably large gravity- gradient drifts 

2 ’ 
from the f /h terms in an inclined orbit. However, in this case, going to the "four ! 

2 '■ 

square" configuration will not help, since then we pick up the terms which are j 

cilso of order 10 mas/yr. In both cases, these terms could probably be modeled and 
removed in the data reduction. Again, ihe polar orbit would avoid this complicated 
proced ure . 


VII. MAGNETIC TORQUES 


The rotating gyro in the gyro experiment will produce the London moment field 
wliich will then interact with its surroundings, including the readout loop, the elec- 
trodes. and any superconducting shields present. This interaction will produce 
back reaction torques which have been analyzed in great detail primarily by Holdeman 
[6.7,8] and Ebner and Sung [9]. A brief discussion and summary of these calcula- 
tions is given here. 

The basic equations governing superconductivity are the London equations [10]. 
Those equations predict that magnetic field is excluded from a nonrotating supercon- 
ductor. except for a very thin region near the surface. Hence, as a good approxima- 
tion, one can assume that the aj^ropriate boundary conditions are that the normal 
component of the magnetic field B is zero at the surface of the superconductor. The 
magnetostatic form of Maxwell's equations in a vacuum is 


A ^ 0 V • A ^ 0 


so that one can define a potential function 1 such that 




^ 0 


2 ! 


and this,, together with the above-mentioned boundary conditions on iS, constitute a 
boundary value problem very similar to traditional electrostatic boundary value prob- 
lems. Solving the problem then allows one to calculate torque on the gyro from the 
integral of the expression for the Lorentz force on the surface currents in the rotor. 

Using this type of approximation the following magnetic torques have been 
analyzed. 

1) Torques due to external magnetic fields on the gyro. These torques have 
been analyzed both with the full London equations and with the approximation men- 
tioned above. The results are similar and are understandable as the traditional 
expression for the torque on a magnetic dipole in a magnetic field where the magnetic 
dipole moment is given by the London moment. 

2) Torques due to the readout ring. These torques have been analyzed with 
and witlvjut a spherical shield. The result is understandable as the torque due to 
the interaction of two magnetic dipoles, one with the London moment and one with the 
magnetic moment of the readout loop. 

3) Torque due to flux coupling. The orientation of the magnetic moment 
effects the flux linked by the readout loop and thus the energy of the system. This 
energy is angle dependent and thus leads to a torque. 

4) Torques due to residual trapped flux in the shield or the electrodes. 

5) Torques due to a flat section in an otherwise spherical shield. 

6) Torques due to holes in a spherical shield. 

7) Dissipation torques froiu electrodes and shields. 

All these torques are found to be negligible provided (1) magnetic field, 
currents, and alignments are kept to levels which are already required for other 
reasons in the experiment: and (2) roll averaging is relied on. These calculations 
rely primarily on simple magnetostatics, so one can have a great deal of confidence 
in them. 

Holdeman [6] has given a simple example to show how a superconducting shield 
can torque the gyro. Consider a plane shield a distance d from the gyro. Then, 
one can use the method of images familiar from electrostatics to insure that the 
boundary conditions are satisfied at the shield. The torque will then be given by 




3 n (itr-n) - 
(2d)'^ 




where 6 is the field of the image magnetic dipole M' and n is a unit vector connecting 
the image dipole with the real dipole, both with magnetic moment equtd to the x^ondon 
moment. For d ■ 4 cm and . “ 200 Hz, we find drift rates of order 
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''' ^ 

where = 1 tnas/yr and a is the angle between n and M*. Tlw superconducting 

shields, by excluding magnetic fields, can produce relatively large drifts, and 
symmetry must be relied on (in this case a 0 or a ir/ 2). 

While a spherical shield will produce no torque on the gyro, a cylindrical shield 
will, and one can get a general idea of the order of magnitude from the plane calcu- 
lation. So if cylindrical shields are used, care must be taken to preserve the sym- 
metry of the situation. This means the gyro must be on the axis of the cylinder and 
either parallel or perpendicualr to it. Otherwise the torques could be unacceptably 
large. The calculation of torque due to a cylindrical shield has not been carried out 
and would require extensive analytical and numerical work. 


VIII. COSMIC RAY IMPACTS 


The energy loss of a minimum ionizing proton per unit length dE/dx is about 
2 

2 meV/gm/cm when traversing material due to electromagnetic interactions. Its 
momentum loss AP is approximately 1/c dE/dx Ax for a relativistic particle where Ax 
is the path length. If all this momentum loss were transferred to the gyro a deflection 
of the spin axis A O would result of order 


2 dK 

JL_ 

I .c" 


- 7.5 ■ 10 rad 


2 

For heavy cosmic rays tins result is multiplied by z where z is the atomic number. 

If N particles were incident on the gyro from random directions, tlien we have a 
random walk process and the total deflection would be multiplied by For protons 


2 

the flux f^^ of particles is about 2/cm /.sec so for At - 1 year 


N ^ C “f At) ^ 6.8 « 10® 
o p 


and 


N 


':1 


- 4 - 10 ® mas/yr 



This is a factor of iO^ smaller than that of Reference 3 Section L3(f) due to a dropping 
of a factor of w in the denominator of the equation in that section. For heavy nuclei, is 

2 

increased by a factor of Z /f /f where f is the flux of particles of atomic number z. 

z p z * 

For iron this can result in about a two order of magnitude increase in Q. For nuclear 


interactions, if n secondaries are produced tl^n is modified by a factor 


‘P 


where f^ is the flux that produces showers and this can result in an increase in Q of 

about an order of magnitude. So we see that cosmic rays are not a problem. Solar 
flares and particles trapped in the radiation belts are also not significant since the 
number of such particles which reach the gyro in a year is not much larger than the 
number of cosmic rays impacting the rotor. Any directional effect would destroy the 


dependence of but this is unlikely and would only increase by about 4 orders 
of magnitude. 


IX. CONCLUSION 


The conclusion was reached that for the current gyro concept of the Gyro 
Relativity Experiment the Newtonian drifts can be reduced to levels well below the 
relativistic effects to be measured provided the conditions on rotor sphericity, mass 
unbalance, and charge and magnetic field, residual gas pressure, and ali^mment 
described in this report are met. These requirements, though difficult to meet, seem 
to be feasible. The most challenging conditions appear to be coating uniformity to 
meet the mass unbalance criteria and residual gas pressure. Also, u sufficiently low 
g environment is crucial for the experiment. Such a low g environment is only avail- 
able in space and is the reason for the dramatic improvement of gyro performance 
possible in the Gyro Relativity Experiment. 
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